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ABSTRACT 

We explored the role of X-ray binaries composed by a black hole and a massive stellar 
companion (BHXs) as sources of kinetic feedback by using hydrodynamical cosmolog¬ 
ical simulations. Following previous results, our BHX model selects low metal-poor 
stars (Z = [0,10 -4 ]) as possible progenitors. The model that better reproduces ob¬ 
servations assumes that a ~ 20% fraction of low-metallicity black holes are in binary 
systems which produce BHXs. These sources are estimated to deposit ~ 10 52 erg of 
kinetic energy per event. With these parameters and in the simulated volume, we find 
that the energy injected by BHXs represents ~ 30% of the total energy released by 
SNII and BHX events at redshift z ~ 7 and then decreases rapidly as baryons get 
chemically enriched. Haloes with virial masses smaller than ~ 10 10 M© (or T v - Ir < 10 5 
K) are the most directly affected ones by BHX feedback. These haloes host galaxies 
with stellar masses in the range 10' — 10 s M©. Our results show that BHX feedback is 
able to keep the interstellar medium warm, without removing a significant gas fraction, 
in agreement with previous analytical calculations. Consequently, the stellar-to-dark 
matter mass ratio is better reproduced at high redshift. Our model also predicts a 
stronger evolution of the number of galaxies as a function of the stellar mass with red¬ 
shift when BHX feedback is considered. These findings support previous claims that 
the BHXs could be an effective source of feedback in early stages of galaxy evolution. 
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1 INTRODUCTION 


The regulation of star formation activity (SF) in galax¬ 
ies of different masses is still an open problem in cos¬ 
mological simulations. Different feedback mechanisms have 
been proposed to regulate the transformation of gas into 
stars. Among them, supernova (SN) feedback is consid¬ 
ered to play a critical role in the regulation of SF, with 
a larger impact for galaxies with circu lar velocities lower 
than ~ lOOkms -1 (iDekel fc Silkl fl986l ) . The SN feedback 
modelling in hydrodynamical codes has certainly helped 
to reach better agreement of the proper ties of the sim¬ 
ulated galaxies with observations (e.g. Governato et al.l 
120071 : ISpringel et al.l 120061 : IScannapieco et ah 20061 ). How- 
ever, these studies also show that SN feedback might not 
be efficient enough to regulate the SF activity in low-mass 
systems, principally at high redshift. This problem mani- 
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fests itself, for example, in the apparent inconsistency be¬ 
tween the stellar-to-virial mass relation obtained from abun¬ 
dance matching techniques and that of hydrodynamica l 
simulations in the low-mass end (e.g. ISawala et al.l 1201 ll : 
iMoster et al.ll2013l : lBehroozi et al.ll2013l ). Different feedback 
mechanisms such as cosmic rays, radiative feedback or pho¬ 
toionization, among others, have been explored with the 
aim of imp roving the regulation of the SF in numerica l 
models (e.g. Wadepuh l fe Springelll2011 - L Stinson et al.ll20ll : 


iHopkins et al.ll2014l : ICeverino et al.ll2014l ). However, much 

work is still needed to both understand the physical pro¬ 
cesses which regulate star formation as a function of redshift 
and to improve their modelling in numerical simulations. 


High-mass X-ray binaries (HMXBs) are accretion- 
powered stellar systems composed by a black hole or a neu¬ 
tron star and a companion massive star, w hich emit X rays 
at typical luminosities of 10 38 ergs _1 (see iFabbiand l200d . 
for a review) and hence provide a radiative feedback on 
the interstellar medium. In some cases, they also produce 
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collimated outflows (jets) with similar kinetic luminosities. 
These jets act as another feedback mechanism, heating the 
surrounding medium. The interest on HMBXs has increased 
in the last years because theoretical and observational works 
suggest that the production, and possibly the luminosity of 
these sources, could increas e with decreasing metallicity of 


the stellar progenitor (e.g. Drav 200d; iLinden et al.l l20ld : 


iKaaret et al.ll201ll : iBrorbv et al]|2014| j. Under this hypoth¬ 
esis, HMBXs could play an important role in the formation 
and evolution of galaxies at early epochs of the Universe. 

Given their radiative feedback, HMXBs have been pro¬ 
posed as effective sources of reionization of the intergalactic 
medium (IGM) and of heating of the gas comp onent in the 
small haloes in the very e arly Universe (e.g. Power et al.l 
l200flh . A recent report from lFragos et ahl (|201 ? ) use the re¬ 
sults of combining semi-analytical m odels of galaxy for ma- 
tion with a population synthesis code dFragos et alj20l3l . to 
compare the X-ray luminosity density produced by HMXBs 
and active galactic nuclei as a function of redshift, claim¬ 
ing that the first prevails at 2 > 6 — 8. Different works 
have also studied the role of these binaries in the heat¬ 
ing of the IG M and the shaping of the 21 cm signal from 
first galaxies JPmyer etal .| [20lli : lKaarelll2()iil : iFialkov et al.l 
120141 : IPacucci et alJl2014h . HMXBs have also been jnvesti- 
gated as a source of reionization. iMirabel et al.l (1201 lh sug¬ 
gest that the rate of ionizing photons emitted by HMXBs 
with black holes might be greater than that of their pro¬ 
genitor stars, and that these sources could heat the IGM up 
to 10 4 K, keeping it ionized. Using zoomed hy drodynami- 
cal simulations of mini-haloes, I Jeon et~aki (j20l4 1 have anal¬ 
ysed the feedback from Population III (hereafter Pop III) 
HMXBs. Their results suggest that X-ray photons from 
these sources may suppress small scale structures and re¬ 
duce the recombination rate in t he IGM, provi d ing a net 
positive feedback on reionization. iKnevitt et al.l (120141 1 in¬ 
vestigated the effect of HMXBs on the high redshift IGM 
using a one-dimensional radiative transfer code. Contrary 
to previous works, they found that HMXBs do not produce 
neither any significant additional ionization nor heating of 
the IGM, except for the distant IGM in the case of con¬ 
tinuous star formation. As it can be seen from the above 
discussion, the issue of the effects of the radiative feedback 
from HMXBs is far from being settled. 

Recently, ljustham fc Schawinskil (l2012h studied high- 
mass X-ray binaries (HMXBs) as a potential source of feed¬ 
back through the kinetic energy of their jets. Those HMXBs 
comprising black holes (hereafter, BHXs) are among the 
most powerful X-rays and jet emitters. The potential of 
BHXs as an efficient feedback mechanis m at hi gh red s hift i s 
supported by observational results from lFender et aTI (l2005l i 
which suggest that the contribution of kinetic energy from 
BHX jets could be significant compared with the e nergy in¬ 
jected by SN, in some cases. I Ramsey et al.1 J2006ll studied 
the interstellar environment of seven HMXBs, showing that 
to understand the ionization and kinematics of the super¬ 
shells around these sources, it is necessary that HMXBs con¬ 
tribute with a significant amount of kinetic energy. Several 
works indicate that the energy deposited by BHXs in the 
form of kinetic energy could be as high as their bolomet- 
ric X-ray lum inosit ies dG allo et al.ll2005l ; IPakull et al.ll20fol : 
iFeng fe Soriall201ll; 

vational evidences, 


Sori^et^I. 20j4|) . Base d on th ese obser- 


Justham fe Schawinskil d2012l 'l conclude 


that the early injection of energy by these sources could heat 
up the ambient gas without expelling it from the galaxy, 
preventing the early transformation of gas into stars, and 
changing the properties of the ISM where SN events will 
take place. As we mentioned before, the main impact of this 
process is expected to occur at high redshift, where low- 
metallicity sources prevail. 

In this paper, we implement for the first time, a BHX 
feedback model within cosmological hydrodynamical simu¬ 
lation, in a self-consistent way. Our simulations in clude a 
chemical evolution model ( Scannap ieco et al. 200514 and a 
physically-motivated SN feedback (IScannapieco et all2006h . 
therefore the enrichment of baryons can be followed as galax¬ 
ies are assembled. Hence, our simulations allows us to de¬ 
scribe the metallicity dependence of the sources along the 
Hubble time. This gives us a further insight into this prob¬ 
lem, by providing the self-consistent evolution of the gas 
cooling, the transformation of gas into stars and the stellar 
evolution, with the subsequent chemical enrichment of the 
baryons. We will consider very low-metallicity stars as BHX 
progenitors, since they are expected to be most energetic 
and abundant. We explore the effect of BHX feedback on 
the cosmic star formation history and on the properties of 
galaxies within haloes of different masses. Our results are 
in agr eement to those reported by ljustham fe Schawinskil 
(120121) where a semi-analytical model is developed to study 
the effects of BHX feedback. 

This paper is summarized as follows. In Sect. [2] we de¬ 
scribe the numerical simulations and the implementation of 
BHX feedback. In Sect.[3]we explore the effects of this feed¬ 
back on the cosmic star formation history, while in Sect. [4] 
we study its effect on simulated galaxies. Finally, we present 
our main conclusions in Sect. [5] 


2 NUMERICAL METHODOLOGY 

In this section, we describe the main characteristics of the 
cosmological code, the simulations, and the model developed 
to incorporate the BHXs self-consistently in the numerical 
code. 


2.1 The cosmological code 

We use an e xtended versio n of the TreePM/SPH code P- 
GADGET-3 dSpringel 120051 4. which includes a multiphase 
model for the gas compone nt, metal-dependen t cooling an d 
SN feedback as described by lScannapieco et al.l ()2005l . l2006l i. 

The multiphase scheme allows the coexistence of low- 
and high-density gas clouds, diminishing the problems of 
over-smoothing of standard SPH. The main difference is 
that neighbouring particles are selected not only by a dis¬ 
tance criterion, but also by considering their relative en¬ 
tropy. The decisions are made on a particle-particle basis, 
which implies no mass-d ependent parameters to be fixed (see 
IScannapieco et al .1120061 . for details). The multiphase model 
works coherently with the SN feedback model, which in¬ 
cludes Type II and Type la SN events (hereafter, SNII and 
SNIa, respectively). When SN events are produced, metals 
and energy are distributed within the hot and cold gaseous 
phases surrounding the star particle representing a single 
stellar population. The cold phase of the star particle is 
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defined by gas particles with T < 2 T, and p > O.lp*, 
while the rest of the gas determines th e hot phase (p t = 
7 x 10 -2 ® gem -3 and T, = 4 x 10 4 K, see lScannapieco et al.l 
120061 for a detailed discussion on these parameters). The 
fraction of energy distributed into the cold phase is defined 
by the parameter e c . An exploration of t his parame ter and 
its effects is given by lScannapieco et al.l (|2006l . [20081 1. Here, 
we adopt e c = 0.5 as it has been done in previous works 
which used this SN feedback model (Scarmapieco et ahll2009l : 
iDe Rossi et alJ [20131 : IPedrosa et al. 12014 1. Note that these 
cold and hot gaseous phases are defined for the stellar pop¬ 
ulations where SN events are estimated to be produced, and 
for the only purpose of distributing metals and energy. They 
have no direct relation with the SPH integration itself. Hot 
gas particles thermalize immediately the SN energy they 
receive, whereas cold gas particles build up an energy reser¬ 
voir. The feedback energy is stored in these reservoirs until 
gas particles have enough energy to change their entropy 
in order to match that of their surrounding hot neighbour 
media. When this occurs, the reservoir energy is pumped 
into the internal energy. This SN scheme is able to produce 
powerful, mass-load ed galactic winds as has been shown in 
previous works le.g. IScannapieco et al.ll2008l . l2009l l. 

Stars more massive than 8 Mq are considered as SNII 
progenitors. We estimate the number o f SN II by adopt¬ 
ing the Initial Mass Function (IMF) of Chabrien (|2003j). 
For SNIa, we use the simple model of lMosconi et al] ( 200il l 
where the SNIa rate is estimated by adopting an observa- 
tionally motivated ratio (~ 0.0015). The lifetimes of SNIa 
progenitors are randomly selected within the range [0.1- 
1] Gyr. This model, albeit simple, allows a good represen- 
tation of the chemic al patterns of the stellar populations 
(I.Timenez et al.l 120141 '). The initial chemical composition of 
the gas component is assumed to be primordial (Ah = 0.76, 
Xne = 0.24). The chemical model follows the enrichment 
by twelve isotopes: 4 H, 2 He, 12 C, le O, 24 Mg, 28 Si, 56 Fe, 
14 N, 20 Ne, 32 S, 40 Ca and 62 Zn. For SNIa we adopt the nu- 
cleosyn thesis model of Iwarnoto et al. (1999) and for SNII 
that of IWooslev fc Weaver! ( 1995lf . Chemical elements are 
distributed in a similar fashion as the energy, although 80% 
of the new elements are dumped into the cold phase (i.e. 
20% goes into th e hot phase). The cooling functions are 
metal-dependent dSutherland fc Dopita|[l993 '). 


2.1.1 Numerical experiments 

The simulated volumes represent 14 Mpc comoving-side 
boxes, and are consistent with the cosmology A-CDM with 
Ha = 0.7, = 0.3, fib = 0.04, <js = 0.9, and Ho = 

lOO/i Mpc -1 kms -1 with h = 0.7. Although this is not the 
current favourite set of cosmological parameters, we use it 
because it allow us to compare these simulations with previ¬ 
ous ones. The variation of the cosmological parameters will 
not affect the conclusions of this work, which is related to 
the effects of feedback on the properties of the ISM and the 
regulation of star formation within individual galaxies. For 
this purpose, we compare runs of the same initial conditions 
but with different feedback mechanisms. 

Initially, the simulations contain 230 3 dark matter par¬ 
ticles and 230 3 gas particles (we will refer to these simu¬ 
lations as S230). The dark matter and initial gas particle 
masses are ~ 9.1 x 10 6 Mq and ~ 1.3 x 10® Mq. The adopted 


gravitational softening length is 1.24 h 1 kpc. S230 runs 
were followed down to z = 0, but we focus the analysis on 
z > 4. 

Higher numerical resolution runs of 2 x 320 3 particles 
were also performed to investigate the dependence on mass 
resolution (we will refer to these simulations as S320). The 
dark matter and initial gas mass in this case are ~ 3.1 x 
10® Mq and ~ 4.9 x 10 5 Mq, and the gravitational softening 
length is 0.5/i _1 kpc. These simulations were run down to 
z ~ 7 due to their high computational cost. We adopted 
the same star formation and feedback parameters used in 
S230. As a consequence, the star formation activity starts at 
higher redshift since smaller systems are better resolved. As 
a result, the impact of BHXs occurs also at higher redshift 
in S320 than in S230. Nevertheless, this behaviour does not 
affect our analysis since we are always comparing runs with 
and without BHX feedback with the same resolution to draw 
conclusions. 

We have run the same initial condition with SN feed¬ 
back (S230-SN and S320-SN) and with SN+BHX feedback 
(S230-BHX and S320-BHX) in order to assess the effects 
clearly. Although in this work, we only discuss the most 
successful implementation, we have run the same initial con¬ 
ditions varying the parameters of the BHX models, as ex¬ 
plained in the following section. 


2.2 The stellar feedback model for BH-HMXBs 

We are interested in studying the effects of the kinematic 
feedback of BHXs originated from very low-metallicity pro¬ 
genitors. As these metal-poor progenitors are expected to 
form principally in the Early Universe, their contribution 
should be relevant for the regulation of the star formation in 
the early stages of galaxy evolution. As mentio ned in the In¬ 
trodu ction, we considered previous r e sults from lPower et al.l 
(|2009l l and [justham fe Schawinskii (l2012fl to develop our 
model, which is fully implemented within our SPH code. 

We will consider BHX populations rather than individ¬ 
ual sources due to the numerical resolution of our simula¬ 
tions. We assume that massive stars with metallicities in the 
range Zemp = [0, 10 -4 ) will give origin to extremely metal- 
poor progenitors of BHXs (see below for more details). This 
metallicity range is taken as indicative, since there is no 
clear theoretical or observational value to be used. Never¬ 
theless, theoretic al works predict lar ger formation rates for 
Z < O.OIZq fe.g. iLinden et aDl2010l l. The upper limit rep- 
resent s the characteristic va lue of the oldest Population II 
stars dBelc.zvnski et af]l2004i '). We would like to stress that 
although the upper metallicity limit was selected with this 
criterion, sensitive variations should not affect strongly our 
results. Note also that we are not modelling Pop III stars. 
We assume our sources to be very low metallicity stars but 
not representative of the first stars. Pop III sources are be¬ 
yond the scope of this paper and would demand the treat¬ 
ment of other p hysical mechanism s which are not included 
in our code (e.g. IChen et al.ll20lj f. Our model includes the 
Z = 0 stellar populations since there are not efficient mix¬ 
ing processes within the SPH kernel. Their treatment should 
increase the metallicity floor quickly from Z = 0. 

In order to estimate the number of BHs produced 
in each young stell ar population, we adopt the model of 
iGeorgy et al.l ( 20091 1 for the evolution of massive stars. This 
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model includes stellar rotation and a dependence of the type 
of compact remnant on the metallicity of the progenitor. The 
fraction of massive stars that end their life as BHs increases 
with decreasing metallicity. Therefore, for each young stel¬ 
lar population we esti mate the n umbe r of BHs produced by 
adopting the IMF of IChabrierl d2003T l. We consider stellar 
masses in the range of 30 — 120 M§ and metallicities consis¬ 
tent with the defined Zemp interval. We assume that ~ 20% 
of low-metallicity BHs will end up in binary systems, pro¬ 
ducing a BHX (/bhx” = 0.20). This value is r oughly com¬ 
parable to the 30 % reported by lPower et al.l J2009l l. after 
correcting by the different IMF adopted. 

As we mentioned before, in BHXs the BH accretes mass 
from its companion star releasing a significant amount of 
energy in the form of radiation and outflows. Observational 
evidence from shock-ionized bubbles in the environment of 
these sources (~ 50-300 pc) suggests that the ir mechani¬ 
cal o utput energy would be around ~ 10 52 erg (ISoria et al.l 
I2014T ). On the other hand, many authors have reported that 
the kinetical power from accreting black holes could be as 
high as their Eddington l uminosity limit dPakull fe Mirionil 
20031: [Pakull et al.ll2Q10l ). For example, for Cygnus X-l 


Gallo et al. ( 20051 ) have found that its jet transforms its to¬ 


tal power into kinetic energy, and that the latter could be 
greater than its bolo metric X-ra y lumi nosity. Similar results 
have been found by ISoria et al. J201C ) for the microquasar 
S26 in NGC 7793, and bv lsoria et alJ ( 20141 1 for an accreting 
BH in M83. 

Motivated by these observations, we assume that the 
energy deposited by low-metallicity BHXs in the form of 
kinetic energy is Ebhx = 10 52 erg dMirabel et al.l l201ll : 
iFeng fe Soria 120 111!. Assuming that a similar amount of en¬ 
ergy is radiated in the X-ray band during ~ 3Myr, the 
mean X-ray luminosity per source is Lx,source ~ 10 38 erg s _1 , 
which is typical of HMXBs. We also tested larger and smaller 
energy values. Energies of the order of 7?bhx = 10 53 erg 
sweep the gas in haloes too efficiently while lower values 
(Ebhx ~ 10 51 erg) produce no significant impact on the reg¬ 
ulation of the star formation activity. 

Hence, for each young stellar population with metallic¬ 
ity in the range Zemp = [0,10 -4 ), we estimate the number 
of BHXs as Nbhx = /Ihx x N B h(Z) where N B h (Z) is the 
number of BHs produced according to lGeorgy et al.l (120091 ) 
and /bhlx = 0.20, as we discussed before. Each of these 
Abhx sources releases Ebhx = 10 52 erg into the ISM. We 
assume that this energy will be efficiently thermalized by 
the surrounding gas clouds where the the events took place. 
This hypothesis is supported by observational evidences that 
show that these events are expected to affect t he region sur¬ 
roun ding the stellar progenitors (~ 300 pc, e.g. lPakull et al.l 
l20ldh As a consequence, the released energy is dumped into 
the internal energy of the nearby cold gas phase, contribut¬ 
ing to build up their energy reservoir. 


3 THE COSMIC STAR FORMATION RATE 
DENSITY 

The SN and BHX feedback mechanisms modify the thermo¬ 
dynamical properties of the ISM, thus affecting the transfor¬ 
mation of gas into stars, which is the main effect we analyse 
in this paper. In order to have a global picture of their im¬ 



Z 


Figure 1. The cosmic star formation rate density estimated from 
simulations with SN feedback (S230-SN, blue dashed line) and 
with both SN and BHX feedbacks (S230-BHX, red solid lin e). For 
comparison, the observations compiled bv lBehroozi et ahl l? 2013i ) 
are also included (black filled circles). The two runs share the 
same cosmology, and the SF and SN feedback parameters. The 
cosmic SFR is diminished by BHX feedback at early epochs (z > 
7) through cold gas heating, and boosted at lower redshift as more 
gas is available to feed the star formation activity. 


pact, we estimate the cosmic star formation rate density 
(cSFR) for both simulations, S230-SN and S230-BHX. As 
shown in Fig. [T] the confron tation of the cSFR w ith obser¬ 
vational results compiled bv lBehroozi et ahl ( 2013 ) suggests 
that the treatment of BHX feedback worked in the expected 
way, yielding an improved description of the observations. 
The factor /bhx ~ 0.20 was actually chosen to be able to 
represent this observational relatiorjj. However, this factor 
might vary if the Ebhx is changed accordingly so that the to¬ 
tal amount of released energy per population remains within 
the same values. Otherwise, the effects on the ISM are found 
to be too weak or too strong. The f act th at it agrees roughly 
with that derived by lPower et al.l (120091 ) is encouraging. 

In order to understand the combined effects of the two 
feedback mechanisms, in Fig. [5] (upper panel), we show the 
comoving energy density (e) released by SNII and BHXs, as 
a function of redshift, for both S230-SN and S230-BHX. We 
only considered SNII contributions since they also originate 
from massive stars. SNIa events will also contribute but with 
a larger time-scale and lower rate. 

Both the SNII and the BHX energy contributions grow 
with redshift, due to the increasing cosmic star formation 
activity. However, BHX energy feedback grows at a slower 
rate than SNII one, due to the metallicity dependence of 
the progenitors. As the metallicity of the ISM increases, the 
number of BHX events diminishes, making SNe the dom¬ 
inating feedback mechanism. As can be seen from Fig. [21 
BHX feedback is significant for z > 7. As an example, at 
z > 7 it represents ~ 30% of the total energy (SNII+BHX) 
released into the ISM. For lower redshift, the contribution 


1 For the adopted 1?bhx = 10 52 erg, we tested between 

~ 1% and ~ 20%, finding ~ 0.20 reproduced the observed 

trend (Fig. 0. 
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redshift 


Figure 2. Upper panel: Comoving energy density released by 
different feedback mechanisms as a function of redshift: SN1I in 
S230-SN (magenta solid line), SNI1 in SN230-BHX (green dashed 
line), and BHXs in S230-BHX (blue dot-dashed line). The to¬ 
tal comoving energy density released by all agents in S230-BHXs 
(BHXs + SNe) is also plotted (red dotted line). BHXs release a 
significant amount of energy, of ~30% of the SNIf+BHX feedback 
at z > 7, but its fractional contribution decays for lower redshifts, 
down to ~10% at z ~ 4. Lower panel: Same as above, but taking 
into account only the energy injected into the cold gas phase of 
the ISM. The contribution of BHXs is larger in this case because 
they dump all their energy into the cold phase. 


of BHX feedback progressively decays, reaching less than 
~ 10% of the total at 2 ~ 4 (~ 5% at z ~ 2). The to¬ 
tal e released by SNII and BHX together is comparable to 
that of S230-SN for z > 7, while for lower redshift, that of 
S230-BXH is greater due to the SFR enhancement triggered 
in this simulation, as can be seen in Fig. [0. We will come 
back to this point later on. 

The SNII and BHX energy release should affect the reg¬ 
ulation of the star formation in different ways, depending on 
the energy available to heat up the cold gas clumps where 
stars are born. Fig. [2] (lower panel) displays the energy re¬ 
leased by SNII and BHXs to the cold gas phase, as a function 
of redshift. For S230-BHX, the amount of energy dumped by 


2 The impact of BHX feedback on the global SF activity when the 
numerical resolution is increased is similar, although as explained 
before, S320 runs formed more stars for the same set of SN and SF 
parameters (see Section 2.1.1). By comparing the energy released 
by BHX and SNII events in S320-SN and S320-BHX, we found 
a maximum of ~ 30% contribution from BHX feedback, that 
decays rapidly with redshift as the ISMs are chemically enriched, 
in a similar fashion as reported for S230 in Section 3. The specific 
redshift at which BHX feedback is maximum over SNII feedback 
is a statement which depends on the box-size. Our simulated box 
correspond to a small field region, so a larger volume will be 
required to make a full statement on this aspect. 


BXHs at high redshift is similar to that of SNII. The total 
amount of feedback into the cold gas is then larger in S230- 
BHX than in S230-SN, producing an early heating of the 
gas that diminishes the star formation activity. In fact, this 
can be appreciated from the cSFR, shown in Fig. [T]for both 
simulations. As expected, in S230-BHX the BHX feedback 
diminishes the star formation activity for z > 7 with respect 
to that of S230-SN, because it provides a larger amount 
of energy used directly in cold gas heating. The lower star 
formation rate in S230-BHX also explains the lower SNII 
feedback in this simulation (Fig. [2]|. We note that the star 
formation activity in S230-BHX increases for z < 7 faster 
than in S230-SN, reac hing values comparable to observations 
jBehroozi et al.ll2013l ). Hence, the impact of BHXs on the 
regulation of the star formation is boosted to lower redshift 
by making gas available to star formation at late times. 


4 THE EFFECTS OF BHX FEEDBACK ON 

SIMULATED GALAXIES 

In order to assess the effects of BHX feedback on the regu¬ 
lation of the star formation in galaxies of different masses, 
we first constructed galaxy catalogues from S230-SN and 
S230-BHX. We use the Friends-of-Friends technique to se- 
lect the virialized st ructures, and the subfind algorithm 
dSpringel et al.l 1200 ll ) to identify the substructures within 
the virial radii. We only consider simulated galaxies resolved 
with more than 500 particles. In order to test the robust¬ 
ness of our results against numerical resolution, we run a 
higher numerical resolution initial condition increasing the 
mass resolution by a factor of eight. We acknowledge the fact 
that our simulated volume is small to provide a complete de¬ 
scription of the galaxy populations. However, by comparing 
the same initial condition run with different feedback agents, 
we are able to underpin the effects of these agents within the 
mass range of ~ 10 9 —10 11 Mq. Hence, the analysis should be 
taken on individual galaxy basis, and as a first step toward 
understanding the impact of BHXs as a possible feedback 
mechanism. 

We carry out an analysis of the simulated galaxies from 
z ~ 9 to z ~ 4. For each available snapshot of the simu¬ 
lations, simulated galaxies are defined at the optical radius 
which corresponds to that enclosing 83% of the baryonic 
mass of the selected system. The dark matter halo masses 
are calculated at the virial radius. 


4.1 The stellar-to-dark mass ratio 


First, we study the ratio between the mean stellar mass and 
their virial halo mass of the simulated galaxies as a func¬ 
tion of the latter. The estimated relations for S230-SN and 
S230-BHX are shown in Fig. [3] For comparison, we also 
show the res ults obtained u s ing th e abundance matching 
techniq ue b y BehroozietalJ d2013h at redshift z ~ 7, 5, 4, 
a nd by Moster_et_ak ( 20131 ) at redshift z ~ 4. In the case 
of lMoster et alj ( 20131 ) data, we extrapolate their stellar-to- 
halo mass relation to lower halo masses. We note that the 
simulations tend t o resolve smalle r viria l haloes compared 
to those shown bv lBehroozi et all d2013l ). except for 3 ~ 7 
where there is a better match of the mass range. Therefore, 
we extrapolate the trends to draw conclusions. 
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At z ~ 7, there is an excess of stars in small haloes 
in S230-SN compared to those detected in S230-BHX. This 
excess originates in the very efficient transformation of gas 
into stars at the early stages of evolu tion in S230-SN as 
also reported i n pre vious works (e.g. ISawala et al.l 1201 ll : 
iDe Rossi et al.l 120131') . Although our SN feedback model 
is successful at producing galactic outflows and fountains, 
which regulate the star formation, it is less efficient in the 
first stages of galaxy formation. When the energy released 
by BHXs is added to the SN feedback, then a decrease in 
the star formation activity at high redshift is detected. In 
Fig. 0 we compared the relative importance of both feed¬ 
back mechanisms as a function of redshift. In Fig. [3]we can 
appreciate the effects on individual galaxies. In fact, the 
regulation of the star formation by BHX events is more effi¬ 
cient in low mass haloes, steepening the relation between the 
stellar-to-halo mass ratio and the halo mass of the galaxies. 
In haloes less massive than ~ 10 10 M©, simulated galaxies 
in S230-BHX have less stellar mass than those in S230-SN. 
For larger halo masses instead, the simulated galaxies in 
S230-BHX reach larger fractions of stars compared to the 
SN feedback run. This trend is produced because simulated 
galaxies in these haloes have larger amount of gas to feed the 
subsequent star formation process. This can be explained as 
a consequence of the impact of BHXs at early times, when 
the systems were smaller and could be affected by BHX 
feedback. 

We note that although the cSFR of S230-BHX repro¬ 
duces the observations compiled bv iBehroozi et all (l2013f) . 
there seems to be still more stars per dark matter halo at 
z < 7. However, the slope of stellar-to-halo mass relation 
agrees better to the extrapolation of observat ions for lower 
masse s. This is also valid when compared to iMoster et al.l 
(120131) at 2 - 4. 

4.2 Gas fractions 

The impact of the SN and BHX feedback on the stellar mass 
of galaxies can be also studied by comparing the ratio be¬ 
tween the mean gas mass and the virial halo mass. Fig. [4] 
shows this relation as a function of the halo mass for galaxies 
in S230-BHX and S230-SN, and for the same redshift range 
displayed in Fig. [3j The gas mass of the simulated galaxies 
is estimated as the total mass of gas particles inside their 
optical radii. 

As can be seen from Fig. [4] simulated galaxies in S230- 
BHX have higher gas masses than those in the simulation 
without BHX feedback, regardless of their host haloes. At 
2 ~ 7, we detect the smaller gas fraction in S230-SN for 
Mh > 40 10 M©. This result seems to be at odds with the 
trends shown in Fig. under the hypothesis of a close- 
box model. Indeed, if there were no gas infall or outflow, 
then the galaxies with the larger stellar mass fraction should 
have had the smaller gas fraction. However, this is not what 
we see in Fig. [4] The fact that galaxies in S230-BHX have 
larger fraction of gas within their optical radii suggests that 
their counterparts in S230-SN have lost larger fractions of 
their gas reservoirs and hence, the subsequent star forma¬ 
tion is lower. If this were the case, then galactic outflows 
should have carried out more material to the circumgalactic 
medium (CGM) of these galaxies, increasing its tempera¬ 
ture. 





[M e ] 

Figure 3. Mean stellar-to-halo mass ratio as a function of halo 
mass, for redshifts z ~ 7, 5 and 4. We estimate the stellar mass 
of a galaxy as the total mass of stellar particles inside the optical 
radius. Dashed lines represent galaxies in S230-SN, while solid 
lines show those in S23 0-BHX. The results of abundance matching 
techniques taken from IBehroozi et alj (120131) are plotted in black 
as open squares plus a dott ed line. For z ~4we include also the 
extrapolation of the data of IMoster et alj (120131) (dash-dotted red 
line). 
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Figure 4. Gas-to-halo mass ratio as a function of halo mass 
for simulated galaxies in S230-BHX (solid lines) and S230-SN 
(dashed lines). The relations are shown for three different red- 
shifts: 2 = 7 (green triangles), 5 (blue squares), and 4 (violet 
circles). Gas mass is estimated as the total mass of gas particles 
enclosed within the optical radius of a galaxy at a given redshift. 


In order to analyse this feedback loop between the ISM 
and the CGM of the simulated galaxies, we estimate the 
fraction of the total gas mass within the virial radius which 
remains within the optical radius in both simulations. As can 
be seen in Fig. [5] in S230-BHX, simulated galaxies in haloes 
with virial masses in the range 10 10 Mg < Mh < 10 11 Mg 
are able to retain larger amounts of gas than their counter¬ 
parts in S230-SN. For smaller and larger haloes, both simu¬ 
lations show the same level of remnant gas within the optical 
radius. In the quoted halo-mass range, galaxies in S230-SN 
have been able to loose a larger fraction of gas, which can 
be understood considering their higher early star-formation 
activity. As a consequence, gas outflows are triggered, ex¬ 
pelling larger amount of gas. This produces a decrease of 
the star formation activity afterwards since there is less gas 
to fuel it. In fact, as shown in Fig. [l] the star formation 
activity of S2 30-SN at low redshift remains lower than that 
estimated bv lBehroozi et aTI d2013l ). whereas it is higher at 
high redshift. Hence, this suggests that it will not be pos¬ 
sible to reach an agreement if only SN feedback is present 
(see also IStinson et al.ll2012l ) . Increasing the SN feedback to 
decrease the star formation at high redshift would make the 
discrepancy at lower redshift even larger. It is the action of 
the BHX feedback that helps to regulate the star formation 
at very high redshifts, which then results in the desired trend 
at lower ones. 

As we can see from Fig. [3] there seems to be a halo mass 
threshold below which BHXs are efficient at modulating the 
star formation. This is approximately Mh ~ 10 10 Mg, which 
corresponds to circular velocities I4ir ~ dOkms^ 1 or a virial 
temperature of T v i r ~ 50, 000 K (assuming primordial abun¬ 
dances). For galaxies below this threshold, the BHX feed¬ 
back is strong enough to heat the ISM to T ~ T V1I , decreas¬ 
ing the star formation activity without expelling a significant 
fraction of gas mass. 



M h [Mo] 

Figure 5. Fraction of the total gas within the virial radius re¬ 
tained within the optical radius, for simulated galaxies in S230- 
BHX (solid lines) and in S230-SN (dashed lines) at z = 7 (green 
triangles), 5 (blue squares), and 4 (violet circles). 


To illustrate this, in Fig. [6] we show the mean tempera¬ 
ture of the gas within the optical (left panels) and the virial 
(right panels) radius, for simulated galaxies in both runs, at 
z ^ 7 and z ~ 5. The run with BHXs shows hotter environ¬ 
ment within the optical radius, although with temperatures 
smaller than T v ; r . Conversely, galaxies in S230-SN have lost 
gas by SN feedback. In fact, SN feedback has been success¬ 
ful at building a hot CGM even at these high redshifts. The 
CGMs are colder in the run with S230-BHX, since the total 
stellar mass formed per halo is not enough to drive power¬ 
ful winds yet. The BHX feedback has contributed to heat up 
the gas within the galaxies, decreasing the mass of new-born 
stars in haloes with Mh < 10 10 Mg. Therefore, as they grow 
by hierarchical clustering, they will be able to have larger 
gas reservoirs to continue their star formation activity. Later 
on, as the SN rate increases, the SN feedback contributes to 
heat up the CGM of galaxies in S230-BHX. 


4.3 The stellar mass of galaxies 

To assess how the effects of BHX feedback modifies the stel¬ 
lar mass of galaxies, we estimate the cumulative number of 
galaxies as a function of their stellar mass. As it can be seen 
from Fig. [7] the simulations behave differently. The simula¬ 
tion with BHX feedback shows a notable evolution with red¬ 
shift, with ~ 80% of the galaxies exhibiting stellar masses 
smaller than 10 7 Mg at z ~ 7, percentage which decreases 
to ~ 50% at 2 ~ 4. Conversely, the simulation with only SN 
feedback does not present a clear evolution in this redshift 
range. On average, ~ 50 — 60% of the galaxies have stellar 
masses smaller than ~ 10 7 Mg. These trends confirm that 
galaxies in S230-BHX form less stars in small galaxies than 
those in S230-SN. Hence, BHX thermal feedback affects di¬ 
rectly small galaxies but indirectly larger ones. Those mostly 
affected have stellar masses between 10 7 and 10 8 Mg and 
inhabit haloes with V™ ~ 40 km s _1 . 
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Figure 6. Median temperatures of the gas component within the optical radius (left panels) and within the virial radius (right panels), 
for galaxies in S230-BHX and S230-SN at z = 7 and z = 5. The dotted lines represent the relation between virial mass and virial 
temperature. Note that at z = 7, simulated galaxies in S230-BHX have hotter interstellar media than their counterparts in S230-SN. 
Conversely, the circumgalactic media of galaxies in S230-SN get hotter sooner, as a consequence of the hot material transported by 
outflows. These outflows are triggered by the more violent SN feedback, produced by the larger SFR taken place in this run at very high 
redshift. 
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Figure 7. Cumulative distribution of the stellar mass within the optical radius at z ~ 7 (green dashed line), z ~ 5 (blue dotted line) 
and z ~ 4 (violet solid line). Left: Simulation without BHX feedback. Right: Simulation with BHX feedback. 
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Figure 8. Mean stellar-to-virial mass relation at z = 7 and z = 
8, for simulated galaxies in S320-SN (dashed lines) and S320- 
BHX (solid li nes). The res ults of abundance matching technique 
from iBehroozi et alj j2013T ) are included (dotted line plus empty 
squares). The main impact of BHX feedback occurs for A/), < 
10 10 Mg, in agreement with results found for S230-BHX. 


4.4 Testing numerical resolution 

In order to test the robustness of our trends against numeri¬ 
cal resolution, we analyse the simulated galaxies in S320-SN 
and S320-BHX at t ~ 7.5 and z ~ 8. Galaxies have been 
selected and analysed by using the same criteria, although 
the minimum number of particles was increased to 2000 in 
order to avoid including very small systems which were not 
considered in S230-BHX. 

In Fig. [5] we show the stellar-to-virial mass relation in 
a similar fashion as in Fig. [3] The relations for galaxies in 
S320-BHX and S320-SN cross each other at about Mh ~ 
10 10 Mg in agreement with the trends found in the S230 
runs. Even more, the cumulative distribution of stellar mass 
also exhibits a similar trend to those found in S230 (Fig. 0. 
And as a consequence, this agreement between the low and 
high resolution runs suggest that the results are numerically 
robust. 


4.5 Contribution of X-ray photons 


As already discussed, the rate and X-ray luminosity of BHXs 
are expected to be highly related to the metallicity of the 
progenitor stars as suggested by observational and theo¬ 
retical studies. Due to the long-mean free path of X-ray 
photons, BHXs could contribute significantly to the heating 
and ioniz ation of the IGM at early stages of galaxy evolu¬ 
tion (e.g. iMirabel et al.l |201l|; |jeon et al.l 20141; IPower et all 


120131 : iKnevitt et al.l 20141) . Mirabel et al. ( 2011 ) show how 
the evolution of the temperature of the low-density neutral 
IGM due to heating by X-rays from BHXs may depend on 
the value of the parameter /x. This parameter is defined as 
the ratio of the total X-ray luminosity Lx of a galaxy in the 
2 10 keV band to its SFR, and its value can be estimated 
from the expression 


/x = (3.5 x 10 40 SFR) -1 L x , (1) 

where Lx is measured in ergs -1 and the SFR in Mg yr -1 . 
An increase of /x would cause the neutral IGM to be heated 
earlier due to the presence of BHXs. This could limit the cold 
gas accretion in dwarf galaxies at high redshift. 

Although the detailed analysis of X-ray feedback of 
BHXs into the IGM is out of the scope of this work, we es¬ 
timate the value of /x using our model for extremely metal- 
poor BHXs in order to evaluate its compatibility with ex¬ 
pectations from other models. The emission of these BHXs 
should dominate the X-ray luminosity of star-forming galax¬ 
ies mainly at high redshift when most young stellar popula¬ 
tions were metal-poor. In any case, our predicted /x repre¬ 
sents a lower limit, resulting in a conservative contribution 
to radiative feedback in the early Universe. 

Assuming that BHXs radiate a similar amount of en¬ 
ergy in X-ray luminosity to that injected in kinematic out¬ 
flows (i.e. ex = -Ebhx = 10 52 erg), and a mean emission 
time of ~ 3 Myr, the mean X-ray luminosity per source as 
Lx ~ 10 38 ergs -1 , as already mentioned in Section 2.2. If 
the number of BHXs produced in a galaxy is IVbhx, then 
the mean total X-ray luminosity of a simulated galaxy orig¬ 
inated from BHXs is Lx = -/VbhxLx, which can be easily 
computed at each time step of the simulation . To make a 
proper comparison with iMirabel et al.l (1201 lh . we use the 
galaxy catalogue of S320-BHX at z ~ 9. These simulated 
galaxies have SFR and L® in the range [0.01 — 20] Mg yr -1 
and [10 39 — 10 41 ]ergs -1 , respectively. Hence, the mean /x 
value obtained at z ~ 9 is (/x) = 0.65. 

Note that the total X-ray luminosities of simulated 
galaxies do not exceed the X-ray luminosity limit of the 
deep est surveys (4 M s Chandra Deep Field-South —CDF- 
S—. IXue et al.l 1201 3) . At z ~ 9, the 2-10 keV rest frame 
energy band corresponds roughly to the observed soft band 
0.2-1 keV. The flux limit of the 4 Ms CDF-S in the soft band 
is 9.1 x 10 -18 ergcm -2 s , therefore the (k-corrected) lumi¬ 
nosity limit at 2 ~ 9 is Lx,ii m ~ 4.6x 10 42 ergs -1 . Hence, our 
BHXs model does not contradict the upper limits imposed 
by current surveys. 

The (fx) value obtained does not take into account the 
absorption of X-ray photons before reaching the IGM. A 
rough estimation of the fraction of escaping photons into 
the IGM could be obtained by adopting a power law spectral 
energy distribution for the sources, with an index of 1.7 in 
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the 0.5-10 keV range dSwartz et al.lliool 'l. We calculate this 
fraction for a typical halo of Mh = 10 s Mq, assuming an ex¬ 
ponential surface density for the galaxy. We assume that the 
ratio of the baryonic-to-dark-matter mass is ma = 0.17, the 
spin parameter is A ~ 0.03 and jd/rrid = 1, where jd is the 
ratio of the angular momentum of baryons to that of the halo 
llMo et al.ll 19981 : iFernandez fc S hulll 1 201 il l . Considering that 
X-ray photons would ionize both h ydrogen and helium in th e 
IGM, we use the cross sections of lKuhlen fe MadarJ (l2005lj . 
The estimated escape fraction is then / eac ~ 0.53. Therefore, 
the effective (/x) would be (/x)eff ~ 0.34 at redshift z ~ 9. 
This value is lower than that proposed by iMirabel et al.l 
J20111) where f eBC ~ 1.0 is assumed. Our simple estimation 
implies that the generated X-ray photons could barely heat 
the temperature of the IGM above 10 3 K. However, our es¬ 
timations are too crude to draw more robust conclusions on 
this aspect. This effec t can be studied with radiative transfer 
models in detail (e.s. iKnevitt et ahlhoMh . 

5 CONCLUSIONS 

In this work, we explore the effects of HMXBs com- 
posed b y an accreting blac k hole on galaxy evolution 
J.Tustham fc Schawinskil l2012h . Following previous results, 
our model assumes that the rate, X-ray luminosity and out- 
flow energy of BHXs increase fo r low metallicity progenitors 
dDravll20f)fil : iLinden et alll20ld l. The BHX feedback model 
is grafted into a version of P-GAD GET-3 which includes SN 
feedback and chemical evolution dScannapieco et ahll200di . 
so that its effects can be studied as the structure forms and 
evolves. 

The estimated kinetic energy released by BHXs is as¬ 
sumed to be efficiently thermalized and hence, is pumped 
into the cold gaseous phase in the surrounding region of the 
sources. This extra source of energy has an impact on the 
properties of the ISM at early stages of evolution and in low- 
velocity haloe s. Encouragingly, our results agree with those 
reported by djustham fe Schawinskil 120121 1 where a semi- 
analytical model was used to implement a scheme based on 
similar hypotheses. 

Our results can be summarized as follows: 

• Following observations which indicate that the kinetic 
energy deposited by these sources might be comparable to 
that radiated in the X-ray bands (during ~ 3 Myr), we esti¬ 
mate that each BHX event could inject ~ 10 52 erg in form 
of kinetic energy in the ISM. With these hypotheses and in 
order to reproduce the observed cSFR, our BHX model re¬ 
quires that a fraction of ~ 20% of BHs with metallicity in 
the range Zemp = [0,10~ 4 ] should end up as X-ray binary 
systems. We tested combinations of larger and lower released 
kinetic energy and BHX fractions finding that they produce 
too strong or too weak effects, resulting in an inadequate 
regulation of the star formation at high redshift. 

• Our model predicts BHX feedback to affect more 
strongly galaxies with stellar masses in the range 10' 
10 s Mq which inhabit haloes with virial velocities smaller 
than V v ir ~ 40kms _1 . Larger haloes would experience a 
negligible impact due to their larger potential wells but 
they could be indirectly affected by the accretion of smaller 
haloes. These results are shown to be robust against numer¬ 
ical resolution. 


• The energy injected by BHXs helps to regulate the star 
formation activity in low-mass haloes, decreasing the frac¬ 
tion of stars formed at very high redshift. As a consequence, 
there is more gas available for star formation at later times. 
Therefore, when BHX feedback is i ncluded, the cSFR re pro- 
duces closely observational results dBehroozi et aljf2013h . In 
this simulated volume, we found the BHX energy represents 
a ~ 30% of total feedback energy released into the ISM at 
z ~ 7, and a 10% at 2 ~ 4. 

• When BHX feedback is included the number of galax¬ 
ies as a function of stellar masses shows a larger evolution 
with redshift (for z > 4) compared to the run with only SN 
feedback turned on. The variation of the evolution is driven 
by the low-mass haloes which are more significantly affected 
by BHX feedback, and hence, its characteristics could be an 
observational test to probe the action of this feedback. 

• The stellar-to-virial mass relation for z ~ 7 is in bet- 
ter agreement to that predicted from abu ndance matching 
dMoster et al.l 120131 : iBehroozi et aid l2013l l. For lower red¬ 
shift, there is still an excess of stars, although the slope of 
the relation is more comparable to the predicted trends. 

• Simulated galaxies with BHX+SN feedbacks exhibit 
hot ISM at very high redshift but their mean temperature 
is smaller than the virial temperature of the haloes. If the 
BHX feedback is switched off, the resulting ISM are colder 
and the fraction of new born stars higher in low mass haloes. 
As a consequence the number of SNe is enough to blow part 
of the ISMs, contributing to building up hotter CGMs. In 
order to prevent this, the SN feedback should be weaken but 
this would lead to a further increase of the star formation, 
producing a larger stellar fraction. 

Hence, our results support previous claims that the 
BHXs could be an important source of feedback in early 
stages of galaxy evolution by regulating the star formation 
in haloes with I4ir ^ 40kms _1 . It goes in the same di¬ 
rection as the so-called ’enhanced’ feedback although BHX 
feedback would have significant impact only in the very early 
Universe. 
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